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Abstract 

We determine the conditions of formation of spontaneous polarization phase transition 
to the superradiance regime in the inverted system of nonrelativistic electrons on equidis- 
tant Landau levels in rarefied magnetized plasma. The possibility of realization of such 
conditions in the lower Jupiter magnetosphere is shown. The effect of cyclotron super- 
radiance on the Landau levels gives a key to interpretation of the nature of superpower 
radioemission of the Jupiter-Io system. 
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1 Introduction. 



The phenomenon of superradiance (SR) was considered for the first time in the known paper 
by Dicke [1] in the example of a two-level model. At present, a significant number of publications 
are devoted to its study (see for example, reviews [2,3,4]), but, as noted by many authors, the 
theme is far from being exhausted, and many interesting and physically important questions 
and situations still remain to be investigated. One of them is the question of the possibility of 
realization of SR in systems with equidistant levels in the presence of inversion. The system of 
fast electrons in a homogeneous magnetic field, for which the transversal movement (rotation) 
energy spectrum is described by the known Landau levels 

E ± =nhcu, n = 0,l,2,..., (1) 

u = — , 2 
mc 

is the most easily realizable and important of such systems [5]. The question on the possibility 
and conditions of SR formation in such systems, besides general physical interest, is also of big 
astrophysical interest because the inversion on the Landau levels is easily achieved in rarefied 
magnetized plasma of active space objects in the presence of bunches of accelerated electrons 
if their initial velocities are directed at some angle to the magnetic field, and so, besides the 
longitudinal energy E», they also possess significant transversal energy E±. For example, in the 
Jupiter magnetosphere [6,7], the observed superpower sporadic nonthermal decameter radiation 
of the Jupiter- Io system with brightness temperature up to 10 15 — 10 17 K can serve as an evidence 
of the possible realization of such situation and generation of SR. 

In the present work, we investigate the question about the possibility and conditions of 
SR formation in the inverted system of electrons on high Landau levels (1) with n ^> 1 and 
the possibility of realization of these conditions in the lower Jupiter magnetosphere. As it 
is known, for the generation of the induced coherent radiation in systems like masers, the 
equidistancy of energy levels is an obstacle because of the specific competition of radiation 
and absorption processes in this case. In the case of SR, we deal not with induced but with 
spontaneous radiation, and here, as we shall see, equidistancy of energy levels appears to 
be of advantage. This is, first of all, because the SR regime is usually realized in open finite 
systems without mirrors, where radiation leaves the active volume of generation quickly enough, 
having practically no time to get in the absorption regime [3]. Second, in this case, all the 
inverted electrons occupy, as a rule, not one level, but some significant interval of high levels 
An (n ^> An ^> 1), and because of equidistancy, all of them radiate the same mode on 
frequency (2), and also, as we shall see, at the same rate independent of the initial energy. 

The phenomenon of SR arises when in "coherence domains", with sizes R smaller then 
the wavelength A, all N Q radiating dipoles gradually during radiation become aligned in one 
direction due to the dipole-dipole interaction between them in the "near zone" (Rq <C A), 
so that, as a result, the total dipole of the domain D turns out to be iVo times larger then 
the elementary dipole a . Therefore, the intensity of the collective dipole radiation becomes 
proportional to Nq, and not iV as in the case of radiation of uncorrelated dipoles. This is 
described in Sec. 2. 

The transition to such a correlated polarized state is similar to the phase transition in 
magnetics or ferroelectrics, and for its description it is convenient to use the Weiss method 
of mean self-consistent field [8]. This theory is developed in Sec. 3. We note that the phase 
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transition under consideration is a nonequilibrium one, and it has all the features of the self- 
organization phenomena in dissipative systems. In Sec. 4, we discuss the key role of the SR 
effect for the interpretation of the observable power of the decameter radiation of the Jupiter-Io 
system. 



2 Cyclotron superradiance on Landau levels. 

For large n, levels (1) correspond to quantum states with wave functions localized near classical 
Larmor orbits with radii Tl = V±/u [5]: 



r . = J4 (2n+1) = ^( 1+ l) R) ^ = ^ = ri . (3) 
V ell V ™>u V 2nJ V rnw 2 uj 

In this connection, we can proceed to (quasi) classical description of such states and transitions 
between them. In a classical limit, in a coordinate system in which the longitudinal movement 
is absent, these orbits are given by 

~f*±(t) = r L {cos(ut + a), sm(u)t + a), 0} , (4) 

where the cartesian components of the vector ~r*j_ are written in braces. 

Being initially inverted on high levels, electrons start to fall down step by step on the ladder 
of states (1), radiating quanta with frequency (2). The differential (angular) and integral 
intensities of the dipole radiation of one electron in the classical limit are described by the 
known formulas [9] 

dl [ d x Tt k ] k 

dn = 4tt C 3 ' Uk = T ] (5) 

7 2eVg 4e^^ 
3c 3 3mc 3 

Equation (6) implies the following evolution law of the electron energy: 

dE±(t) _ _E± _ 3mc 3 

E ± (t) = E ± (0) exp(-t/r). (8) 

We see that the radiation time r does not depend on E±, i.e., electrons with different initial 
energies [within the limits of dispersion AE±(0) E±(0)] will fall down at the same rate. It 
is easy to see that the dispersion of energy will also decrease at the same rate, 

AE ± (t) = AE ± {0) exp(-t/r), (9) 

so that the following relation is satisfied: 

AE ± {t)/E ± {t) = const < 1. (10) 

This allows one to judge about the evolution of the whole collective of electrons by the evolution 
of their average energy and other quantities. Therefore, for simplicity, in what follows, by the 
symbols E±, V±, and d = er L = eV±/u we denote the corresponding values averaged over the 
ensemble of inverted electrons. 
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One can partition the total volume V , occupied by N inverted electrons, into subvolumes 
Vcoh, or "coherence areas" describing the elementary dipoles the sizes i?o of which are smaller 
than the radiated wavelength A but are bigger than the radii of the Larmor orbits (3): 



r L < Rq < A. 

Thus, a large enough number of N radiating dipoles will be in one volume V co h- 

n e V = N^>N = n e V coh > 1. 
Let us consider the total dipole moment of such a subsystem 

N _ N 

%) = £^) = e£r^(t). 

3=1 3=1 



(11) 



(12) 



(13) 



Similarly to (5) for one electron, the collective dipole radiation of our subsystem will be 
described by the formula 



dl = 



x nj] 2 dtt _ u/ 

C 3 47T C 3 



47T ' 



(14) 



By substituting (13) and (4) into this formula and averaging over the period T = 2n/u, we 
have 

(cOs(cjt + CCj) COS^UJt + QLj)) = COs(«j — CXj)/2, 

(sin(cjt + ctj) sin(a;t + ctj)) = cos(o;j — oij)/2, 
(cos(wt + aj) sm(u;t + ctj)) + (sin(cjt + c^) cos(ut + a,)) = 0. (15) 

Furthermore, 

^\t)\ = £ (^(t) • ^(o) = £ (^(0 • ^(t)) + 



No(No-l) _ 

+ £ (^(0.^(0) = 



2 2 

e r r 



JVo(JVo-l) 

A^o + £ cos(«i - a,-) 



and 



(0 ■ Jtf) = e 2 r 
Using the relations 



JVo(JVo-l) 

iV + cos(o;j — ctj) 



(n 2 x + n 2 y )/2. 



(16) 



(17) 



n; 



+ n 2 +n 2 z = l : 1 



+ n J _ 1 + n 2 z 



we get 



(dl) 



e 2 u 2 V 2 



No(No-l) 

N + £ cos(«i - aj) 

i+3 



i + n 2 z dn 

2 4tt' 



(18) 



The factor (1 + ra^)/2 reflects the anisotropy properties of the dipole radiation of nonrelativistic 
electrons. The second term in the square brackets describes the correlation effects connected 
with the mutual aligning of dipoles. If correlations are absent, i.e., if the contributions from all 
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cos(ttj — aj) are mutually compensated and give zero in the sum, then only the first term in 
(18) works, which corresponds to the total radiation of N independent elementary dipoles. It 
is easy to describe situations where the second term in (18) completely compensates the first 
one so that dipole radiation does not arise. These are, of course, the situations where the total 
dipole moment D becomes equal to zero. Consider, for example, the case where N /2 dipoles 
are oriented precisely in one direction, and the rest Nq/2 dipoles are oriented in the opposite 
direction. Then the correlations of pairs in each of these groups are constructive and give 
cos(Aa) = 1; there are N (N Q /2 — l)/2 such pairs in one group and as many in the other, i.e., 
their common positive contribution is equal to N (N /2 — 1). Correlations between pairs formed 
by dypoles of different groups are destructive and give cos7r = — 1 so that the contribution of 
such pairs is, obviously, equal to — 2(N /2) 2 . Adding these contributions together with the first 
term in (18), we obtain: N + N (N /2 - 1) - 2(iV /2) 2 = 0, as it should be from the physical 
viewpoint because we consider the case with D = 0. 

In the case of total correlation, where all cos(a:j — aj) = 1, formula (18) gives 

m _ = N ^Yi l±A £ , (19) 



i.e., the intensity grows in N times as compared with the radiation of N uncorrelated dipoles. 
This is precisely the SR effect [1]. 

The radiation time of such correlated dipoles radiating coherently will decrease No times as 
compared with time (7): 

Tcoh = r/N . (20) 

In reality, it is possible to expect only partial positive correlation of phases, i.e., partial 
aligning of all dipoles, for wich the average over ensemble value of cosine is positive 

(cos(o;j — aj)) = (cos Aa) > 0. (21) 

Having replaced all cos(o;j — aj) in (18) by their average values, we obtain 

p 2 t,) 2 V 2 1 + « 2 <70 

(dl) = [N Q + No(N - 1) (cos Ac)] ™ (22) 

Cr 2 47T 

In this case, the intensity of coherent radiation is proportional to Nq (cos Aa). 

Now we proceed to the consideration of the mechanism of spontaneous aligning of the dipoles 
giving rise to the SR regime. 



3 Polarization phase transition in "coherence 

domains" . 

To solve the problem of the phase transition, we apply here the Weiss method of self- 
consistent mean field [8] confirmed in the theory of spontaneous magnetization. Consider the 
potential energy of a trial dipole "cfo (~^o, t) with the electric field ~$ ("f^o, t) induced at a point 
~r*o by the rest (N — 1) dipoles: 

U(d ) = -~2 ( To, t) ■ ^(Yo, t), (23) 
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where 

E(T ,t) = E U -i 13 — , 24 

All dipoles are rotating according to law (4) and radiate, and their electric field is not static 
but also is rotating with frequency u. Therefore, the use of expression (24) for the field ~fi(t) 
requires explanation. The point is that conditions (11) that determine the "coherence volume" 
imply that different dipoles of a domain are in the so-called near zone (Rq <C A) with respect 
to each other, where the main term in the decomposition of the retarded potentials and fields 
in powers of the small parameters {r L / Rq) and ij- L /\) turns out to be precisely expression (24) 
(in this respect, see, e.g., [9]). 

Averaging Eq. (23) over the rotation period, we notice that the points ~r*o and ~r*j charac- 
terize not the instantaneous positions of the rotating electrons but the positions of the static 
centers of rotation and, consequently, do not depend on time. Substiting (24) into (23) and 
averaging over the period, we obtain 

(U(ro)) = -| E 1 I ^ cos(a - a,-), (25) 



r o 



where do = erL, rij Z = (zo — Zj) /\~7*o — ~f* j\, and the relation S[nj x + n 2 y ] —2 = 1 — 3rij Z is used. 
One can approximate the sum in (25) by the integral with respect to the coordinates rj over 
the "coherence volume" V co h with the obvious measure n e dVj representing the mean number of 
dipoles in the volume element dVj = (d~r*j) in the neighborhood of the point V,. But, prior to 
writing this integral, we notice that, because the position of our trial dipole d (ro,t) can be 
arbitrary, it is necessary to average (25) over this parameter, i.e., to introduce the additional 
integration (dro)/V co h- Moreover, in the spirit of the mean-field method, we replace cos(ao — otj) 
in (25) by its value (cosAa) averaged over the ensemble. After all this averaging, we get the 
expression 

d 2 r r 1 — 1(n ■ ) 2 

(U) = n e (cos Aa) / (dr o )/V coh / dt 3 Lg^ . ( 26 ) 
Z J J r n — T j | 



By making the change of variables {^o^j} {~? = — R = ("^o + "^j)/2} and 
integrating over we get 



d 2 r r 2 - 3^ 2 

(U) = -| n e (cos Aa) J -^(dT>). (27) 

V co h 

Aligning of the dipoles in the same direction is energetically favourable because of reduction 
of the negative contribution to potential energy (U). Therefore, the energetically preferable 
correlations will occur only in that part of the "coherence volume" in which the region of 
integration over the relative coordinates satisfies the condition 

r 2 - 3z 2 > 0, (28) 

i.e., in the region similar to a flattened circular cylinder. We call this part of the coherence 
region by "coherence domain", or by "domain of self-polarization". In a similar neighboring 
domain, the direction of the average vector of polarization should be close to the opposite one 
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to minimize the positive energy of the electric field of polarization in the system as a whole. It 
is known that, by similar reasoning, the macroscopic volumes of magnetics and ferroelectrics 
are also divided into domains. Adjacent domains are separeted by transition regions ("domain 
walls" ) within the limits of which the turning of the polarization vector takes place. 

We return now to the estimation of the integral in (27) with constraint (28). It is convenient 
to calculate it in the cylindrical coordinates (p,z,<p), in which 

r 2 = p 2 + z 2 , r 2 - 3z 2 = p 2 - 2z 2 > 0. (29) 

We consider this integral separately: 

p z - 2z z 



I(pi,p 2 )= Jdtp Jpdp J 2dz { P 2+ * /2 - (30) 



(p 2 + z 2 

pi \r , 

Integration over z in the specified limits gives 2/3V3p 2 and, as a result, 

47T J dp 2 



pi 




7( """^ws/t^ ln| (31) 



(V) = -^=\n{^\ -^McosAa). (32) 
Now we consider the question about the minimal and maximal limits (pi,p2) of the relative 



coordinate p = \J{x\ — x 2 ) 2 + (yi — 1/2) 2 i n the plane between two dipoles in the coherence 
domain. We remind, that the characteristic sizes of the initial "coherence volume" were deter- 
mined by conditions (11): r L <^ R <^ A. Inequality (29) means that, in relative coordinates, 
"the coherence domain" has the form of a flattened circular cylinder with radius 2R . Hence, 
the maximal value of p is p 2 = 2R , and the minimal value of p\ should be taken to be about 
twice the Larmor radius: p\ ~ 2r^, because, at smaller distances between dipoles centers, the 
interaction between a pair of electrons is not of dipole character any more, and it is impossible to 
use the dipole formulas. It must be noted that, under the conditions considered here, the radius 
of the Debay screening is a little bit smaller than r L . Thus, we can write ln(p 2 /pi) ~ ln(-Ro/^i)- 
Because the ratio Ro/r^ enters under a sign of logarithm, the result is weakly sensitive to the 
exact value of this ratio. Thus, taking into account also the condition R <ti A, we can replace 
R here by the quantity of the order of A/10. As a result, we can write 

, / An 2 mc 2 \ , / mc 2 \ . . 

=4w^r ln l^i)- (33) 

Hence, at this stage, inequalities (11) turn into the condition 

E ± < mc 2 /5 w 100 keV. (34) 
With the account of (33), expression (32) takes the following form 

{u) ~ " 371 ln [wj d2 ° Ue (cos Aa) ' (35) 
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To find (cos Aq), we turn to the Weiss method [8]. For this purpose, in the beginning it is 
necessary to consider the response of our system of rotating dipoles in the (x, y) plane 

dj = d {cos(ut + aij), sin(cjt + aij), 0} (36) 

to the external homogeneous electric field 

E? e = E e {cos(ut + cto), sin(cut + a ), 0} (37) 

rotating synchronously with dipoles. The potential energy of a dipole dj in this field, averaged 
over the period of rotation, is 



U = -dj (t) ■ E e (t) = -d E e cos(a j - a ). (38) 

Thus, one can see that the aligning of dipoles along the field with the radiation of released 
energy is energetically favourable. Thermal fluctuations suppress this tendency. They occur in 
rarefied magnetized plasma mainly in the form of plasma fluctuations and Alfven waves. The 
distribution over the phase differences Aa = aj — a is thus given by the Boltzmann formula 

(A ^ r ( ~ U ( Aa ) \ r ( doEeCosAa \ 
p(Aa) = Cexpl — l=Cexpl — I (39) 



with the normalization factor 

2tt 



C- 1 = J d(Aa) exp cos Aa J = 2/ , (40) 

where Iq(x) is the modified Bessel function of zero order [10]. The mean value of (cos Aa) is 
determined by the integral 

2w 2w 

(cosAa) = J cos Aap(Aa)d(Aa) = —In J d(Aa) exp{x(cos Aa)} 
o o 

= ^ h(x) = >m >j!fl (41) 

dx Iq[x) I (x) 

where 

x = d E e /kT. (42) 

So, the external field E e induces the nonzero correlator (41), i.e., in other terms, polarizes 
the system. The measure of polarization is the mean dipole moment of unit volume 

P = n e d (cos Aa) . (43) 

Polarization generates the additional internal electric field 

E p = v ■ P , (44) 

where v is some dimensionless parameter which will be defined below. The field E p , in turn, 
strengthens the polarization. This feedback effect will be taken into account if in (42) we replace 
E e by the sum E e + E p = E e + vP in correspondence with the ideology of the self-consistent 
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mean field developed by Weiss. As a result, we obtain the following nonlinear equation for the 
determination of the polarization P: 

P = n e d (cos Aa) = n e d F ^o(^+ vP) ^j ^ (45) 

where 

F(x)=h(x)/I (x). (46) 
By excluding the external field (E e — > 0), we obtain "the self-consistency equation" for P: 

P = n e d F (^p) . (47) 

We consider now the conditions of existence of its nontrivial solutions. Introducing the 
variable z = {Pd^v /kT), we rewrite equation (47) in the following form: 

* = m- (48) 

The function F(z) has the asymptotics [10] 

F ' 2 » = 1 -i-8F-'"' 2>>1 ' < 49 » 



1-4 + 



z « 1. (50) 



Fiom (48) and (49) it is follows that, at large z, the solution exists and corresponds to the 
polarization of saturation 

* = ^>1. (51) 
kT kT 

Pmax = n e d = — z > — . (52) 
isdo isd 

To determine the threshold value of the electron density above which there arises a nontrivial 
solution of equations (47), (48), it is necessary to consider asymptotic (50). Restricting ourselves 
to the first term, we obtain from (48) 

Z = 2kT Z ' Z ^ ' 

It follows from this equation that the critical (threshold) value of the density n e is determined 
by the condition 

(54) 

For n e > n ec , nontrivial domain self-polarization appears. 

To determine the factor v, it is necessary to compare the expression for potential energy 
following from (43), (44), and (38), 

(U) = —isn e d (cos Aa) , (55) 
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with the previously obtained expression (35). Requiring the equality between them, we obtain 

27r , / mc 2 \ / mc 2 \ . . 



3^ \5E ± 

As an estimate, for example, taking E± = 1 keV, we obtain v m 5.56. 
It is useful to express through the energy E± and magnetic field H: 

2 2 feV ± \ 2 f mcV ± \ 2 2mc 2 E ± 

do = (er L ) = {—) = {-^) = (57) 

As a result, the criterion of the occurrence of domain self-polarization of the inverted electron 
system on high Landau levels leading to SR takes the form 

2tt / mc 2 \ mc 2 n e E\ , . 

^ 1TET )-E7T^>l- (58) 



3V3 \5E ± I H 2 kT 



4 Problem of the power of the decameter radiation of 
Jupiter. 

As an illustration of the application of criterion (58), we consider the values of parameters 
in (58) characteristic for the lower magnetosphere of the active system Jupiter-Io at the base 
of the so-called "Io flux tube", where, according to the observations, one of the sources of 
powerful decameter radiation is located [7]. The characteristic values of the parameters in this 
case are: H ~ 10 Gs, T ~ 10 3 K, and E± ~ 1 keV. Substituting these values into (58), we find 
the critical density of inverted electrons (n e ) c , above which the criterion (58) will be satisfied: 
n e > {n e ) c — 2 • 10 3 cm -3 . This is quite a modest requirement to the density of inverted fast 
electrons at the base of the Io flux tube which can easily be satisfied, so there are all grounds 
to believe that the sporadic super-power decameter radiation of the Jupiter-Io system can be 
connected with the generation of SR. 

Let us consider this question in more detail. We notice at once that, as in the SR regime, 
the energy E±(t) entering (58) quickly decreases, and the density n e satisfying criterion (58) at 
E± = Ej_(0) can cease to satisfy it with the decrease in E±(t). Therefore, for an effective and 
sufficiently long operating SR regime, one obviously requires a large enough excess of n e over 
n ec ~ 2 • 10 3 cm~ 3 . Thus, we conclude that, at the transition to SR regime, electrons at the 
base of the Io flux tube near Jupiter's magnetic pole must have density n e > 10 4 cm -3 . 

According to the basic assumption [6,7] now accepted, the streams of energetic electrons ap- 
pear near Jupiter magnetic poles due to their acceleration up to several keV in the Io ionosphere 
and further movement to the Jupiter along the lines of the magnetic dipole fields. Acceleration 
of electrons is accounted by the induction of electromotive force ~ 400 kV in the Io body and 
ionosphere due to the Io motion through the Jupiter's magnetic field. The lines of the magnetic 
field of the Io flux tube converge to Jupiter's magnetic poles. So, the initial area of acceleration 
of electrons in the Io ionosphere Sq ~ 10 16 -10 17 cm 2 will be "projected" by the magnetic lines to 
an area smaller approximately by two orders of magnitude, Si ~ lO^-So ~ 10 14 -10 15 cm 2 near 
Jupiter's poles. Therefore, the supercritical density of inverted fast electrons n e > 10 4 cm -3 , 
necessary for the occurrence of SR, will be provided by the initial density n > 10 2 cm" 3 near 
Io. It is quite an acceptable value in view of the known data about Io and its ionosphere [11]. 
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Electrons moving along the Io flux tube initially have density below the critical value and, 
therefore, radiate rather weakly But, near the Jupiter magnetic poles, the density of electrons 
reachs the supercritical value, and they pass to the SR regime. Thus, in this region, the 
intensity of collective cyclotron radiation grows approximately by ten orders of magnitude 
(N = n e ■ V co h ~ 10 10 ), what accounts for the observable huge power of decametric radiation 
from a rather small volume of the circumpolar area occupied by the source [7]. A separate work 
will be devoted to more detailed calculation of this power and also to the interpretation of the 
rather unusual dynamical spectra of the decameter radiation of the Jupiter-Io system within 
the framework of the present model. 
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